Gyroscopic effects in interference of matter waves 
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Abstract 

A new gyroscopic interference effect stemming from the Galilean translational factor in the 
matter wave function is pointed out. In contrast to the well-known Sagnac effect that stems from 
the geometric phase and leads to a shift of interference fringes, this effect causes slanting of the 
fringes. We illustrate it by calculations for two split cigar-shaped Bose-Einstein condensates under 
the conditions of a recent experiment, see Y. Shin et al, Phys. Rev. Lett. 92, 050405 (2004). 
Importantly, the measurement of slanting obviates the need of a third reference cloud. 
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^From the spinning mass gyroscope to modern light and matter wave interferometers 
is a spectacular evolution in man's technology for sensing a change in orientation. From 
navigation of ships and spacecrafts to studying the seismic motion of tectonic plates, the 
applications of gyro keep motivating further scientific and engineering progress in the field. 
It thus merits to point out in this Communication a new gyroscopic interference effect that 
stems from the Galilean translational factor in the matter wave function. 

The rotation-sensitive interference experiments have a long history. A rot at ion- induced 
difference in optical paths for two light beams counter-propagating around a closed contour 
causes a shift of the interference fringes, which was conceived by Sagnac as a subject of 
interferometry [jj. This first demonstration of the Sagnac effect was performed on a rotating 
table. Soon afterward, the first measurement of the Earth's rotation was done by Michelson 
a,d Gale 3 with a very ,ar g e station hght interferometer operatin g on the same principle. 
The advent of lasers ensued a rapid progress in the optical measurements of the Sagnac effect 
and culminated with the creation of miniature ring laser gyros currently in use for navigation, 
see a review jj). The Sagnac effect for matter waves was theoretically predicted in The 
result can be obtained by changing light to matter via the substitution fko = mc 2 into the 
expression for the Sagnac phase shift. The advantage of matter waves is in the possibility of 
considerably down-sizing the interferometer owing to the reduced wavelength. The effect was 
demonstrated experimentally first with neutrons 0, 0| and later with neutral atoms 7] and 
charged particles (electrons) S . There is also a suggestion for hybridizing light and matter 
waves by way of polaritons |2|. Highest precision and stability to date were achieved with 
atom interferometers: on a rotating table, at rates of rotation of the order of Earth's rate 
Q, and with stationary apparatus, detectin g the rotation of Earth Q. It can hardly be 
overemphasized that Bose-Einstein condensates (BECs) provide a new type of matter wave 
sources. Their high brightness and coherence make them to matter wave interferometry 
what lasers are to optical interferometry. A single measurement of interference with BECs 
should suffice to determine the rotation velocity; this paper explores BEC as a gyroscope 
for this specific reason apart from a further miniturization. The first demonstration of 



interference between two BECs was reported in 1121. Interference of BECs in lattices was 
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also observed [13, |14J. Recently, two cigar-shaped BECs in a double-well potential were 
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made to interfere keeping the relative phase of atomic clouds under control [15[. The authors 
conclude: "Propagating the separated condensates along a microfabricated wave guide prior 
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to phase readout would create an atom interferometer with an enclosed area, and hence with 
rotation sensitivity" . Thus rotation-sensitive BEC interferometers should become available 
soon. 

In all the above experiments, done and planned, rotation reveals itself via a shift of 
interference fringes owing to the Sagnac effect. In this paper we show that there is another 
effect that leads to slanting of the fringes. 

As noted above, the discussion to follow deals with BEC as a bright source of matter 
wave, so that we base our consideration on the Gr oss-Pitaevs k ii (GP) ^nation Q, bnt the 
argumentation applies also to a dilute atomic cloud [that is including the limit of (3 = in 
Eq. (fT|)]. The phenomenon we present here is thus of a rather general nature, but our specific 
reference to BEC as a matter wave source inevitably leads us to include dynamical effects 
of the nonlinear term. The positive internal pressure expressed by the nonlinear term is 
conducive to rapid expansion of the system when the trap is switched off. The GP equation 
in a rotating reference frame reads 



Here f2 is the angular velocity of rotation, V(r, t) is the trap potential, (3 = 47rNh 2 a/m 
is the coupling constant, where m, a, and N are the atomic mass, scattering length, and 
the number of atoms, respectively, and the order parameter ip{ r ,t) is normalized to unity. 
We shall assume that Q is small and does not depend on time. Let us consider a simple 
gedanken experiment to detect rotation, first with the aid of an approximate representation 
of the experimental steps, and then Eq. (P) will be numerically solved subject to realistic 
experimental conditions. The experiment proceeds in four steps, see Fig. ^ (1) preparation 
of the initial state, which we assume to be the ground state of a BEC in a single-well trap, (2) 
coherent splitting of the BEC into two or more components by deforming the potential into 
several separated traps, (3) independent transportation of the components along some paths 
by varying the positions of the traps, and (4) switching off the potential and observation 
of the interference pattern when the different components begin to overlap as a result of 
their free expansion. Assuming that steps 2 and 3 are done adiabatically slowly the effects 
of rotation on the interference pattern can be deduced analytically. The adiabaticity means 
that the system remains in the ground state at each moment; we will return to this point in 
the discussion of numerical results. During step 3 the solution to Eq. (JTJ) can be decomposed 
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FIG. 1: General scheme of a rotation-sensitive interference experiment with BEC. A BEC, initially 
in a single cigar-shaped trap (1), is split coherently in the transverse direction into two parts (2), 
which are then moved in the longitude direction and released to expand freely (3). The interference 
pattern (4) is obtained by measuring the atomic density when the clouds begin to overlap. 



as 



(2) 



where the overlap of the different components ip n (r,t) is geometrically marginal but physi- 
cally non-negligible due to tunneling. The position of the n-th component is given by 

«(r,i)r^nM)>„ 



(3) 



(V*(r,t)V*(r,t)) n 

where (•••)„ denotes integration over the region occupied by the n-th trap. In this region 
V(r,t) = V n (r — a n (t),t), where V n (r,t) defines the shape of the n-th trap. If a n (t) and 
V n (r, t) are slowly varying functions of t, in the adiabatic approximation we obtain 



V> n (r,t) = exp 



ih 1 fj,(t') dt' + i^nit) 



0n(r;t), 



where and 4> n (r;t) constitute the ground state solution of the stationary equation 



f- - 12 • 1 + V n (r - a n (t), t) + (3\4> n (r; t)\ 2 - u(t) 
2m 



0n(r;t) = 



that depends on t as a parameter, and 7 n (t) is the geometric phase |17l ]. 

ln(t) = % 



' mr;t')dMr;t')/dt') n 
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(0;(r;t')^(r;t')), 

Note that the chemical potential [i(t) does not depend on n and corresponds to the ground 
state of the whole system, which assumes that equilibrium is attained at each instant of 



time via tunneling between the components [la] . Our key observation is that the solution 
to Eq. (jSJ) can be presented in the form 

fi(t) = fi(t) + 0(Q 2 ), (7a) 
<p n (v- 1) = exp [ih^mvnit) ■ r] n (r - a n (t); t) + 0(ft 2 ), (7b) 



where fi(t) and 4> n {r;t) constitute the ground state solution of 

*(r;t) = 0, (8) 



^ - ft ■ 1 + K(r, t) + /3|^(r; t) | 2 - /Z(t) 
2m 



and v n (t) = f2 x & n (t). The exponent in Eq. (|7b)l is the Galilean translational factor; it 
accounts for the fact that the n-th component moves with the velocity v n (t) with respect 
to a non-rotating inertial frame whose origin coincides with that of the rotating frame. 
The solution of Eq. (jHJ) may have a nonzero phase 0(Q) which, however, disappears if 
the potential V n (r,t) becomes axially symmetric about the direction of f2. The validity of 
Eqs. (J7|) can be confirmed for the special case of a harmonic oscillator trap. Substituting into 
Eqs. (JHJ) and (jHJ) V n (r,t) = Vno( r ) — |mr T w 2 r and a n (t) = a, where Cj = diag (uj x , u y , u z ) is 
a diagonal matrix, r is treated as a column, and r T denotes the corresponding row, one can 
obtain an exact relation between the solutions to these equations, 
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/i = /2 — + VHo(b-a), (9a) 

0(r) = exp [ifi _1 mu • r] <p{r — b), (9b) 

where u = Q x b and b = (oj 2 + QQ T — Q 2 ) 1 lj 2 &. It can be easily seen that b = a+0(T2 2 ), 
so the result © is in agreement with Eqs. ((Zj). ^From Eqs. (jSJ)? ©, and (ffbj) we find 

j n (t) = h~ l m / fl x a • da, (10) 



where the integral is calculated along the path traced by the vector a n (t). Summarizing, we 
can specify Eq. © as follows 

-t 



ip(r, t) = exp 



-ih l J p(t')dt' exp [i ln (t)+itT l rrw n {t) • r] n (r- a„(t);t). (11) 



This function provides an initial state for free expansion during step 4. The resulting inter- 
ference pattern depends on the phases of the components in Eq. ([lljh and hence is sensitive 
to rotation. A difference between the geometric phases leads to a shift of the interference 



fringes; this is the Sagnac effect. The Galilean factors introduce an additional phase differ- 
ence; this phase is inhomogeneous, which causes a deformation of the interference pattern. 
There is also a third effect: the interference pattern rotates with respect to the rotating 
frame, similarly to the rotation of the plane of motion of the Foucault pendulum. It is 
the last two effects that may result in slanting of the fringes, as is shown below, and their 
manifestation should be clearly distinguished from the Sagnac effect. 

We illustrate this general consideration by calculations for a two-dimensional system 
modeling the experiment discussed above. Number as well as phase fluctuations in BECs are 
supposedly problematic in measuring a shift of the fringes, but do not affect their slanting. So 
we use the GP equation without explicit representation of the fluctuations 3|. The results 
reported below were obtained from Eq. ((TJ directly, with no approximation or ansatz to the 
GP equation, by propagating the solution in time in a finite spatial rectangle —d < x < d, 
—1/2 < y < 1/2, with the condition ip(x,y,t) = on its boundaries. A BEC containing 
iV = 10 4 atoms 23 Na (m = 3.82 x 1CT 26 kg, a = 2.75 x jttm) is initially in a single-well 
trap, see Fig. [H The trap is formed by a Gaussian potential Vo(x) = —Uq exp(— ^mu 2 x 2 /Uq), 
in the x direction, and zero boundary conditions at y = ±1/2, in the y direction, with the 
trap depth Uq = h x 5 kHz, frequency to = 2tx x 615 Hz, and length I = 23 /xm. In the first 
stage of the calculations, the ground state in the trap is obtained by propagating Eq. (JTJ) 
in imaginary time. The atomic density in this state has a cigar-like shape with the width 
6.8 /iiii on the level of 1% of the maximum value. In the second stage, the atomic cloud is 
split in the x direction into two equal cigar-shaped parts by varying the potential according 
to V(x,t) = f(t)[Vo(x — vt) + Vq(x + vt)] with a constant velocity v during time t x = 5 
ms, where the factor f(t) is introduced to preserve trap's depth at each moment equal to 
Uq, with /(0) = 1/2 and f(r x ) ~ /(oo) = 1. The final distance between the cigars is 
d = 2vt x = 13 /xm, and their width is 5.6 /im, so to a very good approximation (i) cigars 
do not overlap, and (ii) the zero boundary conditions at x = ±d do not affect the solution 
of Eq. Similar coherent splitting of a cigar-shaped BEC into two parts was realized in 



a recent experiment 15 1. The above values of to, Uq, d, and r x coincide with the values of 

these parameters in 15j, but in order to facilitate the calculations the cigar length I, and 

lence the number of atoms N, are taken to be equal to about one tenth of their values in 
~l 

151 ] . In the third stage, the cigars are transported in the y direction with the same constant 
velocity v during time r y . This is done by introducing into the Hamiltonian in Eq. (0) an 
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FIG. 2: An interference pattern obtained by simulation of the experiment shown in Fig. Sunder 
the conditions similar to that realized recently without rotation (a) and with rotation (b). In 
the latter case the interference fringes are shifted (the Sagnac effect) and slanted; the white dashed 
line is a fit to the central fringe according to Eq. (|15|). 

additional term —vp y . The transportation distance s = vr y is chosen to be equal to cigar's 
length, s = 23 //m, which corresponds to r y = 17.7 ms. Finally, in the fourth stage, the 
potential is switched off and the atomic clouds are allowed to expand freely during time 
t = 1 ms. The boundaries at x — ±d still do not affect the results because this time is too 
short for the waves reflected from them could reach the central region \x\ < 4 /im, where 
the interference pattern is observed. The calculations were done without rotation, SI = 0, 
and with rotation for Q = 10~ 2 x u; the results are shown in Fig. 121 Because of finiteness 
of the velocity v , some oscillations of the atomic density are excited during the splitting 
and transportation processes that distort the interference pattern. Indeed, because of the 
cigar shape, the density oscillation along the longitudinal axis has a low frequency and is 
susceptible to nonadiabatic excitations. 

We have repeated the calculations for 20 times slower motion (r x = 100 ms, r y = 354 
ms) with the same values of all the other parameters, see Fig. 01 This situation is much 
closer to the adiabatic regime, and the interference pattern in this case is perfectly smooth. 
As can be seen from the figures, the question of adiabaticity is of considerable importance. 
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FIG. 3: (a) and (b): similar to Fig. [21 but for 20 times slower motion during steps 2 and 3 in 
Fig- E (c): same as in (b), but skipping step 3; there is no Sagnac shift in this case, but slanting 
remains unchanged. 

Elementary excitations caused by nonadiabatic processes would quite generally complicate 
the interference pattern and make the identification of fringes difficult. Hence the experiment 
must be conducted slow enough so that energy imparted to the BEC should not exceed the 
excitation threshold. Anyway, in both cases one can clearly see that in the presence of 
rotation the interference fringes are shifted and slanted. We note that in real experiment 
the phase shift may suffer from uncontrollable offset from run to run while the slanting 
angle is known a priori to be zero in the absence of rotation. An interesting possibility 
of observing the effect of slanting is thus assured by the very manifestation of interference 
fringes. 



To interpret these results, following 



151 ] we assume that the expansion can be described 



by Eq. with (3 = 0. In the adiabatic regime, the initial condition for this equation has 
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the form ([lip. Then for the atomic density in the interference region we obtain 



n(x, y, t) = N\ip(x, y, t) | 2 oc 1 + cos — x' - j s - — y' 




) 



(12) 



where (x, y) and (x', y') are the coordinates of the point of observation in the rotating and 
non-rotating reference frames, respectively, and the time t is measured from the moment of 
release. The first term in the argument of the cosine function in Eq. (|T2|) corresponds to free 
expansion in the x direction of two point-like sources initially located at x — ±d/2. The 
second term comes from the geometric phases in Eq. (fTTj) and is given by 



comes from the Galilean factors in Eq. (fTTj) . Substituting into Eq. (|T2*j) x' = xcos(Qt) — 
ysin(Qt) = x — Qty + 0(Q 2 ) and y' = xsin(Qt) + ycos(Qt) = y + 0(Q), we obtain 



Thus in the presence of rotation the interference fringes are shifted by 73, which is the Sagnac 
effect, and slanted, with the phase shift between their upper and lower ends (y = ±//2) equal 
to 27g- These simple equations work surprisingly well. Indeed, according to Eq. (|TK|) the 
spatial period of the fringes is A = hr/md = 1.35 /mi, which is very close to A calc = 1.38 fim 
obtained by fitting the numerical results. ^From Eqs. (|T3~j) and ()14|) we have 7s = 7g = 4.15 
(recall that in our case s = I). Fitting the position of the central fringe in Fig.|2fb) we obtain 
7g alc = 3.85 and 7q 1c = 5.07; similarly from Fig. 0(b) we find 7g alc = 4.04 and 7 G alc = 4.13. 

To conclude, we have discussed a gyroscopic effect that stems from the Galilean trans- 
lational factor. In contrast to the Sagnac effect that leads to a homogeneous shift of in- 
terference fringes, the Galilean phase causes an inhomogeneous shift between the ends of 
the fringes, i.e., their slanting. For the situation considered (cigar-shaped geometry with 
s ~ /) both effects are of the same order of magnitude. However, to produce a shift two 
steps are needed, transverse separation and longitude transportation, while for slanting the 
first step is sufficient, see Fig. Ofc). Another important difference is that the Sagnac phase 
is subject to uncertainty due to a possible asymmetry between the traps and uncontrollable 



7s = h mVtds. 



(13) 



The third term, where 



7g = h mttdl, 



(14) 




(15) 
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variations from run to run of the experiment which need to be compensated for by using a 
third reference cloud, while slanting is not affected by these factors. 
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